code:

- fit_experimental_data.m:  Fits experimental time-activity data to the biokinetic model to generate the time-activity curves (TACs).

- calc_TCP_metrics.m: Computes the Tumor Control Probabilities (TCP) and metrics of a set of TACs.


data:

- experimental_data.odt: experimental time-activity data for each of the radiopharmaceuticals obtained from different publications.


- fitted_TACs (folder): .mat files that contain the sets of fitted TACs for each radiopharmaceutical (Cu64, Y90, I131) and different time steps (dt = 0.1 h and dt = 0.05 h). The name of the files follows the structure "fitted_isotope_dt.mat". Each .mat file contains a number of Nx2 arrays, which contain N time (column 1, in hours) and relative activity (column 2) points for each curve. There are 7 curves for Cu64, 8 for Y90 and 10 for I131.


- metric_TCP_results (folder): .mat files that contain the metrics and TCP values for each radiopharmaceutical dose rate curve,  sub-lethal damage repair half-life(T1/2) and alpha/beta ratio (3, 10 Gy) combinations. The name of the files follows the structure "isotope_T1/2_alpha/beta.mat". Each .mat file contains:

	· a: scaling factor to convert TACs into dose-rate curves (in Gy/h).
	· alpha, alpha_j: central value and heterogeneous distribution of the parameter alpha of the LQ model (1xN_stoch array, in Gy^-1).
	· beta, beta_j: central value and heterogeneous distribution of the parameter beta of the LQ model (1xN_stoch array, in Gy^-2).
	· G_array: distributions of the protraction factors (MxN_stoch array, where M is the number of curves studied).
	· m1_array, m2_array, m3_array: metrics 1, 2 and 3 values for each TAC (1xM).
	· mu, mu_j: central value and heterogeneous distribution of the sub-lethal damage repair rate (1xN_stoch array, in h^-1).
	· N_cel, N_cel_j: central value and heterogeneous distribution of the total number of cells of the tumor (1xN_stoch array).
	· N_stoch: number of random samples to generate heterogeneous distributions of parameters.
	· rho_1, rho_2, rho_3: Spearman's rank correlation coefficients between TCP and metrics 1, 2 and 3, respectively.
	· sigma: standard deviation considered to generate heterogeneous distributions of parameters.
	· t_rep: sub-lethal damage repair half-time (in h)
	· TCP_array: Tumor Control Probabilities for each curve.


- m1_max_deltaTCP_0.02.mat = metric values that yield a deltaTCP of 0.02 for each radiopharmaceutical, sub-lethal damage repair rate half-life alpha/beta ratio. The file contains 3 sets of arrays, with structures:

	·  "T_rep_array": 1x5 array that contains the values of the sub-lethal damage repair half-times (T1/2 = 0.1, 0.5, 1, 2 and 3 h).
	·  "m1_max_isotope_alpha/beta": 1x5 arrays that contain the values of the metric that yield the maximum deltaTCP established for each of the T1/2 values and alpha/beta value.
	·  "sigma_m1_max_isotope_alpha/beta": 1x5 arrays that contain the uncertainties of the above quantity.

- figure7.mat = data for Figure 7. Contains 9 arrays, 3 for each radiopharmaceutical, with the same quantities explained above, extended to larger T1/2 values (different for each radiopharmaceutical).
